Soft particle model for block copolymers 
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A soft particle model for diblock (AB) copolymer melts is proposed. Each molecule is mapped 
onto two soft spheres built by Gaussian A- and B-monomer distributions. An approximate analytical 
expression for the joint distribution function for the distance between both spheres and their radii 
of gyration is derived which determines the entropic contribution to the intramolecular free energy. 
Adding a mean-field expression for the intermolecular interactions, we obtain the total free energy 
of the system. Based on this free energy, Monte Carlo simulations are carried out to study the 
kinetics of microphase ordering in the bulk and its effect on molecular diffusion. This is followed by 
an analysis of thin films with emphasis on pattern transfer from walls with a periodic structure. It 
is shown that the level of coarse graining in the soft particle model is suitable to describe structural 
and kinetic properties of copolymers on mesoscopic scales. 
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1. Introduction 

Coarse graining by elimination of irrelevant degrees of 
freedom is a central problem in many areas of statisti- 
cal physics, especially when modelling diffusion and the 
phase kinetics in polymer systemsii An important class 
of processes is phase separation in polymer melts with in- 
compatible components, which occur on time and length 
scales many orders of magnitude larger than those related 
to the motion of individual monomers. A detailed un- 
derstanding of such processes is required in present-day 
attempts to utilize polymers in the design of new materi- 
als and in the production of tailored micro- or nanoscale 
structures by self-organizationJ^ 

From the theoretical viewpoint this poses the prob- 
lem to eliminate internal degrees of freedom of polymer 
chains, and to seek for a coarse grained description of the 
kinetics in terms of a small number of collective degrees 
of freedom. In order to cover effects of chain length, 
monomer interactions, heterogeneities in monomer se- 
quences e.g. in the case of copolymers, and interactions 
with walls in confined systems^ such coarse grained the- 
ories should be based on appropriate microscopic input. 

A promising approach is to represent polymer chains 
in their coil conformation in dense systems simply by in- 
terpenetrating "soft particles" . This soft-particle picture 
should be applicable for length scales comparable to or 
larger than the radius of gyration Rg and for time scales 
of the order of the diffusion time td ~ Rq/D, where D 
denotes the center of mass diffusion constant. The soft 
geometrical object assigned to a polymer can be quite 
general and should be characterized by only a few pa- 
rameters. In the simplest case, polymers are described 
by soft spheres with effective interactions that can be de- 
rived by an explicit elimination procedure for monomer 
degrees of freedomi^ Slightly more complex geometrical 
objects are soft ellipsoids, which allow one to take into 
account oricntational effects and anisotropic shape defor- 
mations. For homopolymers the first soft-ellipsoid model 
was developed in Refi^ based on microscopic input de- 



rived from a bead spring model. This model includes the 
excluded volume effect into the effective intramolecular 
free energy. As a consequence, tuning of interaction pa- 
rameters with control parameters (temperature, degree 
of polymerization) is necessary to capture correctly the 
screening effect responsible for the quasi-ideal chain be- 
havior in dense systems. It was shown that this adaption 
of effective interaction parameters can be avoided by us- 
ing the Gaussian chain model as basis for the soft ellip- 
soid model, ^ at the expense that the excluded volume ef- 
fect in this model has to be included in a mean field type 
manner in the self-interacting term of the free energy. 
Based on the Gaussian ellipsoid model, structural and ki- 
netic properties of both one-component melts and binary 
mixtures were studied successfully for bulk systems. In 
a subsequent extension to confined systems it was shown 
how spinodal decomposition in polymer blends becomes 
modified in thin films, including situations of periodically 
patterned walls 

Soft particle models become increasingly valuable in 
the description of complex matter systems, and also ana- 
lytical techniques can be used to relate the input param- 
eters to microscopic propertiesj^iS 

In this paper we propose a soft particle model for block 
copolymers. Phase separation in these systems proceeds 
under the constriction of a chemical link between in- 
compatible polymer blocks FiSiii Depending on the rel- 
ative amounts of different blocks, a variety of ordered 
microphases can emerge on length scales that are tunable 
by the respective degrees of polymerization. Because of 
these unique features block copolymers offer a promis- 
ing tool for the fabrication of nanoscale devices as, for 
example, by the nanolithography technique J^ii^ 

Following the early fundamental work on the deriva- 
tion of a Landau theory by Leibleri* and of fluctuation 
corrections by Ohta and Kawasakii^ and subsequently 
by Fredrickson and Helfandj-- a progressive understand- 
ing of the phase diagram of diblock-copolymers with 
variable length ratio of the two blocks has been devel- 
oped. Important features like the static structure fac- 



2 



tor and the scaling of the microphase periodicity with 
chain length were examined by the lattice Monte Carlo 
technique.— li^ii^ Further simulation studies provided in- 
sight into the detailed changes of chain conformation 
near phase boundaries i^i^i^^ the internal energy and 
entropy,^^ as well as transport properties, e.g. center- 
of-mass diffusion of chains and their modification by the 
order-disorder transition. ^^'^^ More sophisticated meth- 
ods are the "dissipative particle dynamics" (DPD) that 
includes hydrodynamic forces between effective beads of 
polymers,—!^ and methods treating the entangled state 
of copolymers 

On the other hand, for studying the far-from- 
equilibrium dynamics and large scale structures, evo- 
lution equations for monomer densities have been de- 
veloped. One approach is to formulate time-dependent 
Ginzburg-Landau (TDGL) equationsj^i^aiSai^a A pow- 
erful extension thereof is the selfconsistent (SFC) field 
methodf^i^ which proved successful in predicting large 
scale phase morphologies. 

Complementary to these developments is the soft par- 
ticle model proposed here. Motivated by the results of 
Refs.— we explore the possibility to map the internal 
molecular degrees of freedom onto only a few parame- 
ters, allowing straightforward physical interpretation. In 
other words, we attempt to coarse-grain as far as possible 
while keeping the most important structural characteris- 
tics of individual molecules. Because of its simplicity, this 
model should enable one to treat the time evolution of 
ordering structures, including surface-induced structures, 
in an efficient manner. Unlike TDGL or SCF theories, 
it retains orientational effects and shape fluctuations of 
individual molecules, features which again will be impor- 
tant when treating surface effects. Moreover, a model 
of this type could form the basis for treating the phase 
behavior of more complex molecules. 

Our paper is organized as follows. First, we develop in 
section 2 a soft-particle model of (AB)-diblock copoly- 
mers, where each block is represented by a Gaussian 
sphere. Their radii of gyration Rx {X = A, B) and 
the vector r connecting their midpoints are the only pa- 
rameters related to the internal degrees of freedom of a 
molecule. Test calculations for the structure factor, the 
ordering transition, the scaling of the lamellar distance, 
stretching of chains, as well as diffusion properties in the 
bulk are presented in section 3. Thereby we demonstrate 
that this Gaussian disphere model (GDM) reproduces 
known bulk properties in a satisfactory manner. This 
motivates us to treat in section 4 equilibrium and ki- 
netic properties in thin film geometries. Major questions 
concern the wall-induced microphase separation both for 
neutral walls and walls that attract one component, and 
the dynamics of pattern transfer from prepatterned walls 
into the film. Finally, some concluding remarks and di- 
rections for further research are presented in section 5. 



2. Gaussian disphere model for block copolymers 

Following the basic idea in previous work^ we pro- 
pose here to represent each block in an (AB)-diblock 
copolymer molecule by one soft sphere with radius of 
gyration R-^ {X = Aov B), as illustrated in Fig. [T] The 
molecules' orientation is given by the vector r = f*^ — 
connecting the two centers of the spheres. Its magnitude 
r = |r| determines the stretching of the molecule under 
AB-repulsion. The three quantities i?^, R^ and r are the 
parameters that represent the internal degrees of freedom 
of one molecule. In order to express configurational fluc- 
tuations of a molecule in terms of these parameters, we 
have to specify 

i) the joint probability densities P{'r, R^, R^) for 
given number of bonds in the A and B block, Na 
and Nb, respectively. 

ii) the conditional monomer densities pxix\r, R^, R^) 
of the block X written in terms of the coordinate 
X relative to the center of the sphere. The nor- 
malization condition is 

J px{x\f,R^,R'')d^x = Nx + l. (1) 

These densities are to be derived from a microscopic 
model of an isolated chain. For the free energy of an 
ensemble of M molecules with fixed values Na and Nb 
we write 

P — -^intra ~t- -Winter ; (2) 

where 

M 

Fintr. = -kBTY,^nP{n,Rf,Rf) (3) 

i=l 

is the intramolecular part in the absence of intermolec- 
ular interactions. For the latter, contact interactions of 
strength exY among the different blocks {X, Y = A oi 
B) are considered here, 

M 

i^inte. = i E E J2 '^yb' r'yp^ (y) pJ (y) ' 

= l X=A,B Y=A,B •' 

(4) 

where is a contact volume. The parameter b sets 
the microscopic length scale of the model and is used 
as length unit, b — 1. The density pf denotes the X- 
monomer density in the laboratory frame, 

pf (y) = Px{y- r.hn , Rf, Rf) ■ (5) 

Note that we have included the self-interaction terms 
i = j 'm eq. (|4]), which implies that we should consider 
entropic contributions only in Fintra- 

To complete the definition of the free energy, we must 
specify the input functions P and px- Following Refi^, 
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we propose to use Gaussian chains as microscopic input, 
with a few additional approximations that lead to sim- 
ple analytic expressions. First we consider P. Here our 
assumption is that and are uncorrelated, with 
distributions Px{R^) obtained from homogeneous Gaus- 
sian chains. Thus we can write 

P{r,R^,R'')^PA{R^)PB{R'')W{r\R^,R'') (6) 

with 

for large Nx ■ An accurate representation of the scaling 
function is^ 

where K[){z) denotes the modified Bessel function of or- 
der zero and the coefficients a and d are known from 
fits to the exact moments of second and fourth order 
in the limit of large Nx {{{R^?) ^ Nx/6; {{R^)^) ^ 
(19/540)7V|-). For the conditional probabihty in ([5]), we 
make a Gaussian ansatz with respect to f, 

where {r^)ii — 2[(i?'^)^ + (i?^)^] is the variance of j'for 
given i?'^ and R^ . 

As a test of © we compare with Monte Carlo sim- 
ulations the conditional probability W{r\R^) obtained 
from ([5]) by integrating over the orientations of r with 
r = fixed, and by averaging over R^ . Regarding its 
dependence on R^, we took averages over each of eight 
successive intervals on the R^ axis chosen such that they 
have equal weights 1/8 with respect to the distribution 
([7]). Results are shown in Fig.[5]by the continuous curves, 
together with the simulation data. With increasing R^, 
peak positions move to the right, showing that chain re- 
alizations with large R^ imply larger stretchings r. The 
semilogarithmic plot in the inset confirms satisfactory 
agreement between the data and the Gaussian approx- 
imation. 

Second, we turn to the conditional monomer density 
px- As mentioned already, we assume a spherical shape 
and ignore any dependence on r and on the size of the 
opposite block, i.e. px = px{x\R^). Then for large Nx 
and for typical R-^ (excluding highly stretched chains) 
one can show that the scaling form 

,10) 

with f{v) = (3/27r)3/2exp(-3z;2/2) well agrees with 
Monte Carlo simulations for isolated chains. Using (fTU]) 
as approximation in the GDM directly follows the anal- 
ogous reasoning in the soft ellipsoid model, described in 



detail in Ref."S. To complete the description of the GDM, 
we give the expression for the radius of gyration Rq of 
the total chain, to be derived by simple geometric con- 
siderations 

Rl = fA[R^? + /B[i?^]' + fAfsr^ , (11) 

where fx = Nx/N, and N ^ Na + Nb- It is clear that 
this type of model can be extended to more complicated 
structures like triblock copolymers, but we will not pur- 
sue this here. 

So far, for given external parameters Nx, M and V, 
we have constructed a free energy F that depends on 
the variables Rf and r/^ (giving fi = — ff). To 
model the kinetics, these varaiables are changed in a 
Monte Carlo process. Two types of elementary moves 
are considered: (i) block translations rf-^rf+ Arf 
where the components of Ar* are drawn from a uni- 
form distribution in the interval [— Ar^jjx/2, Ar^jj^/2] 
with Ar,^^^ = (3/4)(A'^x/6)^/^, and (ii) size changes 
i?f ^ R^ + Ai?f , where Ai?f is uniformly distributed 
in the interval [-Ai?,-^^J2, Ai?^^j2] with Ai?^^, = 
0.5(7Vj(:/6)"^/^ under the additional constraint Rf > 0. 
All elementary steps have the same a priori probabil- 
ity. Acceptance probabilities are chosen according to the 
Metropolis algorithm based on the free energy 

Results are presented in terms of the Flory-Huggins 
type parameter x — {'2^ab — ^aa — ^bb) / k^T and we will 
use k^T as energy unit in the following, k^T = 1. All 
calculations are performed for caa = £bb = ^ and a total 
mean density of monomers p = 0.85. Unless otherwise 
said, symmetric mixtures with a fraction /a = 0.5 of A 
monomers are considered and a chain length N = 120 is 
chosen. 



3. Bulk properties 

In this section we test the GDM against bulk proper- 
ties of chain-like copolymeres, a subject well known from 
extensive model studies We 
show that despite its simplicity the GDM accounts for a 
remarkable set of both structural and diffusion properties 
connected with microphase separation. In all bulk sim- 
ulations periodic boundary conditions are used and the 
volume V of the cubic simulation box is chosen to con- 
tain in total M = 4000 molecules. Averages are typically 
performed over 10 independent simulation runs. 



3.1 Structure 

A direct measure of short- and long-range ordering is 
the static structure factor S'(fc), which we define through 
the density correlation function of A-monomers, 

S{k) = J d'x {pA{x)pA{Q))e''-' ■ (12) 
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In Fig. [3l data points of the spherical average of ((T2)) are 
displayed for symmetric chains in the disordered phase. 
Continuous curves are fits to the expressioi*^^^ 



with 



Fix) 



NS-\k) = -[FikRa) 
a 




(13) 



(14) 



Equation generalizes the form of the structure fac- 
tor from Lcibler's random phase approximation, which in 
our notation amounts to setting a = 1/2, Rq = Rq and 

5 — 2xN. In this theory, lamellar ordering sets in at the 
critical value (x-^)c = 10.5, connected with a divergence 
of ^(fc) at A; = fc* with fc*i?G = 1-95. Following previous 
.^pj.jji7ji8j22 jj)^ ^ g^j.g regarded as fit parameters, 
which we determine from the behavior of S{k) around its 
main peak at k^,. Their values are given in Table 1. In 
the vicinity of the peak and towards small k a good fit is 
achieved, whereas at larger k the data points fall below 
the continuous curve, see Fig. [3] This is to be expected 
because the intramolecular connectivity is not taken into 
account by assuming Gaussian monomer densities within 
each block. 

The physical meaning of the parameters a, Rq, and 

6 is to allow for a shift in the peak position k^, relative 
to the Leibler value^^ and a deviation of the maximum 
value S{ki,)/N from scaling with x^. Indeed, our results 
in Fig. [3] reflect a downward shift in fc^ with increasing x, 
starting with ki,Rc — 1.8 for x = 0. Qualitatively, this 
downward shift agrees with previous calculationSfi^!^ 
but the effect is smaller than in these works. More- 
over, when plotting NS~^{k^,) from calculations with 
different N against x^i see Fig. |4l linear extrapola- 
tion to zero yields an instability of the disordered phase 
at a critical value (xA^)c that increases with decreasing 
chain length TV. This trend is qualitatively consistent 
with the Fredrickson and Helfand theory^i^ predicting 
(x^)™ ~ {x^)c °^ N~^/^. In our model the critical 
value for the longest chains {N — 300) is roughly esti- 
mated as (xA^)c — 12 to 14. 

Ordered structures spontaneously forming under a 
quench from the disordered state to a value x-^ above 
but still close to {xN)c, display a multidomain pattern 
and weak segregation, i.e. a smooth variation of the re- 
spective monomer densities when passing from A-rich to 
B-rich domains. Figure [S^ exemplifies the structure fac- 
tor for x^ = 30. The position of the sharp peak at 
/c* = 2tt/X reflects the lamellar periodicity A. Obviously, 
ki, is smaller than in the disordered state of molecules 
with the same chain length, due to their stretching under 
alignment. In the GDM, chain stretching is described by 
the parameter r — {f^y^'^. Its dependence on x-^ (with 
N fixed) is plotted in Fig. [SI showing a pronounced in- 
crease in the vicinity of the ordering transition and a sub- 
sequent weaker increase as X-^ rises further. For x = 0, 



one recovers = 2[((i?^)^) + ({R^)^)], see section 2. 
The averaged radii of gyration ((i?"^)^)^/^ of individual 
blocks are practically unaffected by changing x- 

Going to strong segregation, larger scale oriented 
lamellae with sharp interfaces develop. This leads to 
higher-order peaks in S{k) seen in Fig.[S}3, where we plot- 
ted the structure factor for x-^'^ = 54. Note the absence of 
the second order peak due to the arrangement of blocks 
inside the domains. We also carried out bulk simulations 
for asymmetric chains. The existence of a cylindrical 
phase is exemplified by Fig. [5j: for Ja = 0.3, reflected by 
superstructure peaks as marked in the figure. 

Computations with different combinations of x and N 
in the strong segregation range allow us to test the scal- 
ing prediction for the lamellar distance A = 27r/fc* of the 
form X/N^/^ ~ (x-/V)". This relation is obtained by min- 
imizing the sum of the elastic energy due to molecular 
stretching and the interfacial energy within the period A, 
yielding n = 1/6.^ Indeed, a good collapse of data for dif- 
ferent N and x is achieved in the double logarithmic plot 
of \/N^^^ versus x^j as shown in Fig. [71 The exponent 
deduced from the slope of the straight line is n ~ 0.22, 
in fair agreement with the mean field argument indicated 
before. For small (xA'^)-values corresponding to the dis- 
ordered phase, a saturation is seen, i.e. A cx N^^'^. 



3.2 Diffusion 

As shown previously ) "'^^1^^ the center of mass diffusion 
of molecules in the disordered phase slows down upon 
increasing x-^, because of the increasing degree of short 
range order, refiected by the increase of S{k^,). In the 
lamellar phase, diffusion parallel to the interfaces is quite 
similar as in the disordered phase near the ordering tran- 
sition, whereas perpendicular diffusion gets strongly sup- 
pressed and drops to zero when x-^ is increased further. 
As seen in Fig. [8] the GDM reproduces theses general fea- 
tures. To distinguish Dy from D± in the ordered phase, 
oriented lamellae were prepared through appropriate ini- 
tial conditions. The slight increase of Dy right above 
the transition point may originate from the reduction of 
structural fluctuations and smoothening of the interfaces. 
The overall diffusion coefficient is consistent with the av- 
erage D « (2i:»y + D^)/3 « (2/3)L'y. Analysis of the 
time dependent center of mass mean squared displace- 
ments showed that the directionally averaged short-time 
diffusion coefficients are almost unaffected by the inter- 
action X and are close to Z3(x = 0). 

An intuitive approach to understand the diffusion 
across the domain boundaries in the lamellar phase is to 
adopt the picture of one-dimensional Brownian motion 
in a periodic potential V{z). We start from the exact 

expression^ D 



Do/ 



-V{z) „V{z) 



where Dq is the 



diffusion coefficient in the corresponding homogeneous 
state, and the bars denote averages over one period. Ob- 
viously, the rate determining diffusion steps are those 
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near the maximum of V{z) or the minimum of the equi- 
hbrium density Peq{z) oc exp{—V{z)). Within mean field 
theories^ this aspect can be generahzed to interacting 
systems by regarding V{z) as potential of mean force, 
defined in terms of the actual equilibrium density via the 
Boltzmann factor. Guided by these ideas, we write for 
our system 



Do 



(15) 



with simulated A- monomer densities pa{z) and diffusion 
constant Dq at the ordering transition point. Close to the 
ordering transition pa(z) can be represented as pAiz) — 
pj(l + Ap^sinfc^2:), which yields 13^ = Do{l ~ Ap\y^^ . 
This simple approach already describes the sharp drop 
in D± for > {x^)c, shown in Fig. [HI It still over- 
estimates the simulated data for D± especially at larger 
xN. One reason may lie in the fact that shape deforma- 
tions of molecules during barrier crossing are not included 
in these arguments. The barrier crossing problem for a 
block copolymer has been studied in Ref.— . 



4. Thin films 

Confinement of block copolymers between pla- 
nar walls adds several new aspects to microphase 
separation A common repulsion of A and B 

blocks by neutral walls will favor parallel orientation of 
molecules and therefore can induce perpendicular lamel- 
lar ordering. However, if walls act differently on A and 
B-monomers, one type of blocks will be preferred to the 
walls, favoring parallel lamellar ordering. These com- 
peting behaviors, including their time dependence under 
initial conditions corresponding to a disordered state at 
high temperatures, have been studied before. Hence we 
limit ourselves to a few representative and supplementary 
results. 

Patterned walls are particularly interesting for appli- 
cations, where a chemically pre-patterned surface should 
be translated into a correspondingly patterned polymer 
film.^— i^i^ We show that the GDM can provide new in- 
sight into the process of pattern translation. Calculations 
of the time-dependent structure factor enable us to dis- 
cuss the propagation of stripe patterns into the film and 
their competition with spontaneous ordering fluctuations 
away from the pre-patterned substrate. 

In the film simulations periodic boundary conditions 
are used in the lateral directions. Unless otherwise spec- 
ified, films have a lateral size = Ly = 4A with A the 
lamellar distance in the bulk (for example, A — 31 for 
N = 120 and x = 0.45, see section 3). Averages are 
performed over typically 10 independent runs. 



4.1 Homogeneous walls 

We start with planar walls parallel to the {x, ?/)-plane, 
with a potential for A and B monomers as in Refi^, 



yd) _ 



(1) 

ey exp 



^(2) 



exp 



2h 



2/ 



(16) 



where the superscript (1) refers to the lower wall at z = 
and the superscript (2) to the upper wall at z = Lz- 
The quantity 1^, is a parameter that tunes the softness 
of the confinement.^^ For simplicity we have chosen it to 
be the same for A and B monomers and the same for 
both walls, = 0.5. The exponential form of the wall- 
monomer interaction is convenient, since it allows one to 
get analytical expressions for the integrals that determine 
the interaction with Gaussian monomer densities.— 

As discussed above, for neutral walls, molecules near 
the wall acquire a preferential parallel orientation. To 
show this, we set e^^ = e^' = e„ = 1, j = 1,2, and 
calculated the quantity (|cos0(2;)|) in the equilibrium 
state (after thermalization), where Q denotes the angle 
between the vector r of a molecule and the z-axis. The 
results are shown in Fig. [5] after averages were taken over 
X and y for fixed z coordinate. Obviously, angles O near 
tt/2 prevail near z — and z — L^. A preferred paral- 
lel orientation persists even up to the middle of the slab 
(random orientations would give (| cosO(z)|) = 0.5). The 
repulsive interactions between monomers and the walls 
lead to a deformation of molecules close to the walls, 
as seen by the behavior of the average block radius of 
gyration versus z. At this point the neglect of the in- 
ternal molecular structure in the GDM becomes appar- 
ent. Close to a wall, chain polymers have been shown to 
get compressed in the z-direction, but elongated in the 
(x, ?/)-plane.^^ Such anisotropic deformations are beyond 
the scope of the GDM. This sets a lower limit of film 
thicknesses that can reasonably be treated within the 
GDM, i.e. Lz should be larger than a few times the radius 
of gyration. 

To further quantify the order in the film and its evolu- 
tion following a quench from a random initial state, we 
introduce the time-dependent lateral structure factor 



S{k\\,z,t) 



-(|pA(fc||,^,t)P) 



MNa 

defined in terms of lateral density fluctuations 

PA{k\\,z,t) = [ (fx\\pA{x,ty^r-\\ 



(17) 



(18) 



with a; II = {x,y). The prefactor Lz in (|17p is intro- 
duced in order to achieve -independence of (fT7|) in the 
bulk limit — + oo for given total monomer concentra- 
tion. Averaging (fT7|) over with /cy = {k^ + kyY^^ 
fixed and subsequent averaging over z yields the quantity 
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S'(fc|| , t) which is plotted in Fig. [9j3. It clearly reflects lat- 
eral ordering induced by the alignment effect described 
above, with some initial coarsening. During the coars- 
ening the position of the maximum k^^(t) approaches a 
non-zero value, corresponding to the lamellar distance in 
the equilibrium state. The structure is in fact close to the 
structure found in the bulk in Fig. O). In particular, a 
3rd-order peak appears, reflecting a well ordered lamellar 
structure perpendicular to the walls. 

The situation changes when considering an A- 
attractive wall, caused by a stronger repulsion for B than 
for A monomers, 

+ (19) 

where = 1 as before and — —S:^ = S^^ > 0. Choos- 
ing — 0.5, the energetic preference of A-monomers 
is strong enough to overcome the essentially entropic 
alignment effect discussed before that molecules near the 
wall have reduced orientational degrees of freedom. In 
Fig. llOb . the corresponding A- monomer density pa{z) is 
plotted across a fllm of thickness — X (with A from the 
corresponding bulk simulations) . Layers of A-monomers 
adjacent to the walls rapidly form and get separated by a 

B-rich domain. When the wall at 2 = is replaced by a 

(2) (2) 

neutral one (e)^ — ~ 1)' equilibrium density pro- 
file essentially remains symmetric, see Fig.lTOb. Thus the 
A-attractive wall induces layering almost as in Fig. [TUb . 
Development of the A-rich layer near the neutral wall, 
however, takes much longer time than in Fig. llOb . The 
whole situation is reminiscent of wall-induced spinodal 
decomposition in films of binary polymer blends, at least 
in its early stagesj^i^i^ 

4.2 Periodically patterned walls 

In this section we study the kinetics of pattern transfer 

into the film. The wall near 2; = is assumed to be 
(2) (2) 

neutral {e\ = — 1), while the wall near z — has a 
modulation along the y-axis with period Lp, 

e^^\y)^e^[l + 6^cos{kpy)]. (20) 

Here = 2, — —(5^ — Sw — 0.5, and kp = 2'K/Lp. 

As long as X < Xc the equilibrium state displays a peri- 
odic segregation of monomers in the y-direction confined 
near z = 0. The decay of the segregation amplitude along 
the z-axis rcfiects the correlation length of the disordered 
phase. On the other hand, for x > Xc various scenarios of 
pattern-induced microphase separation emerge, depend- 
ing on film thickness and on the commensurability be- 
tween the two length scales Lp and A. 

Let us begin with the commensurate case Lp/X — 1 
and a thickness — X. Figure [TT] shows, for x — 0.45, 
the time evolution of the integrated structure factor 
S{ky,z,t) = J dkxS{k\\, z,t). After a rapid initial growth 
of a peak with ky ~ kp near the patterned substrate, 



the associated periodic structure propagates along the z- 
direction. Simultaneously, a broader structure develops 
across the slab and indicates spontaneous microphase or- 
dering similar to bulk behavior. ''^ At t = 2 x lO'' MCS, 
the wall-induced sharp structure dominates everywhere 
and has penetrated the slab nearly uniformly. A weak 
3rd-order peak ist clearly seen at ky ~ 3fcp. The result 
is an equilibrium state with nearly perfect perpendicu- 
lar lamellar ordering. No changes are observed any more 
when going to larger times. 

For a thicker film with ^ 1.8A, see Fig. [H "bulk- 
like" ordering processes throughout the slab evolve fur- 
ther and more strongly interfere with the ordering wave 
propagating from the patterned wall. For long times the 
ordering wave finally succeeds to overcome the bulk-like 
domain structures. Equilibrium is nearly but not fully 
reached at t = 4 x 10^ MCS. 

Now we turn to incommensurate situations. The case 
Lp/X — 1.1 and — X will be described without show- 
ing explicit results because of its similarity at long times 
to the commensurate case Lp/X = 1.0 shown in Fig. [TT] 
The main difference at shorter times is most directly ob- 
served from the structure factor 5(^11,^) defined as in 
Fig. [5]3. In the commensurate case bulk ordering ini- 
tially leads to a shoulder in S{k^\,t) next to the primary 
peak at kp, which merges with it as time proceeds. On 
the other hand, for Lp/X = 1.1 and i = 200 MCS one 
can identify a separate "bulk" peak. It also merges with 
the peak at kp for longer times t > 2 x 10^ MCS so 
that the final structure corresponds to complete penetra- 
tion of the wall pattern, as in the commensurate case. 
It is interesting to note that transient coexistence of two 
peaks and final dominance of the peak at kp have been 
observed recently in two-dimensional Fourier transforms 
of real space images of diblock copolymer thin films on 
chemically nanopatterned substrates 4^ 

By contrast, for the larger ratio Lp/X — 1.2 with 
Lz = X the two main features in Figs. 113b - c pre- 
vailing near the patterned wall and deeper in the film, 
respectively, remain separated with respect to ky as time 
proceeds. Up to i = 2 x 10'^ MCS the peak at kp domi- 
nates, which is evident from the averaged structure factor 
S'(/fc|| , t) plotted in Fig. [H whereas for t > 2 x 10^ MCS 
the quantity S'(fc|| , t) takes its maximum at the bulk value 
ko = 27r/A. No further change of this pattern is observed 
in our longest runs up to t = 6 x lO'* MCS, nor does it 
change when we increase the strength in the modulation 
of the wall potential to 6^ = 0.75. Hence it appears that 
no well-ordered structure develops in this case, showing 
that transfer of the wall pattern to the film sensitively de- 
pends on the commensurability of the two length scales 
Lr, and A. 



5 Summary and Outlook 

The Gaussian disphere model (GDM) provides a highly 
coarse grained description of diblock copolymer melts at 
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the molecular level which still captures the essential fea- 
tures of self-organized structure formation. The essence 
of the model is to parameterize the internal molecular de- 
grees of freedom by a few stochastic variables, the radii 
of gyration Ra and Rb oi each block and the distance 
vector r . Molecular positions together with these internal 
variables move according to a kinetic MC algorithm. The 
algorithm is driven by a free energy F{{fi}, {Rf}, {Rf}) 
derived from Gaussian chains, in analogy to the Gaussian 
ellipsoid model for polymer melts proposed earlier 

Regarding bulk ordering and diffusion, several features 
known from less coarse grained models were reproduced 
to a good approximation. This prompted us to study 
microphase ordering in thin films, especially the kinetics 
of pattern transfer from a stripe-patterned substrate into 
the film. Detailed results for the time-dependent struc- 
ture factor S{ky, z, t) were presented, including some dis- 
cussion of commensurability effects and pattern penetra- 
tion into films of varying thickness. 

The GDM is expected to loose reliability, when the 



film thickness becomes less than a few times the radius 
of gyration. Then the assumed sphericity of individual 
blocks is no longer compatible with the molecules' dis- 
tortion under confinement. Improvement on that issue 
could be achieved in the spirit of Ref.^ by allowing el- 
lipsoidal block shapes or by representing each block as a 
string of soft spheres. 

On the other hand, the GDM may turn out advanta- 
geous in a semi-quantitative description of more complex 
molecules, including chain-like segmented or branched 
copolymers, or copolymers carrying reactive groups^ 
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TABLE L Parameters for fitting the generalized Leibler func- 
tion (|13p to structure factor data in Fig. [3l 
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FIG. 2: Comparison between the conditional probability VF(r|7?^) based on the Gaussian approximation Q (continuous curves) 
with Monte Carlo data for symmetric chains with A'^ = 100. To display the dependence on Ra, averages have been taken over 
eight successive (7?^)-intervals, see text. Different curves refer to the first five of these intervals. With increasing _R'*-values, 
distributions W(r\R'^) shift towards larger r-values. The inset shows, in a semi-logarithmic representation, the results for 
the smallest and largest of these eight (_R^)-intervals, confirming that the Gaussian approximation is satisfactory as long as 

r <iVl/2. 
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FIG. 3: Simulated structure factor S{k) in the disordered phase for different x^- Continuous curves are fits to Eq. (|13|l . 
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FIG. 4: Normalized inverse maximum of the structure factor versus xN for three different chain lengths. 
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FIG. 5: Spherically averaged structure factor in ordered phases. Lamellar phase with (a) = 30 (weak segregation) and (b) 
xN = 54 (strong segregation), displaying the third order peak, (c) Cylindrical phase with /a = 0.3, = 65 (iV = 100), with 
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FIG. 10: (a) Time evolution ol the A-monomer density in a film of thickness Lz = \ with A-attractive walls, for x = 0.45 (b) 
Same, but with A-attractive left and neutral right wall. 
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FIG. 11: Time evolution of the structure factor S{ky,z,t) in the presence of a stripe-patterned wall near z = for x ~ 0.45. 
The pattern periodicity is Lp = 27r/fcp = A and the film thickness Lz = A. The lateral system size is ~ Ly — 4Lp. (a) 
t = 200 MCS, (b) t = 2 X 10=* MCS, (c) t = 2 x 10* MCS. 
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FIG. 12: Same as Fig. [TTl but L-^ = 1.8A, and t = 4 x 10'' MCS in (c). 
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FIG. 13: Same as Fig. [TT] but Lp = 1.2A, and f = 6 x 10* MCS in (c). 
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FIG. 14: Circularly and z-averaged structure factor from the same simulation data as in Fig. 1131 



